Antiferromagnetic transitions in 'T-like' BiFeO^ 
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Recent studies have reported the existence of an epitaxially-stabilized tetragonal- like ('T-like') 
monoclinic phase in BiFeOs thin-films with high levels of compressive strain. While their structural 
and ferroelectric properties are different than those of rhombohedral-like ('i?-like') films with lower 
levels of strain, little information exists on magnetic properties. Here, we report a detailed neutron 
scattering study of a nearly phase-pure film of T-like BiFeOs- By tracking the temperature depen- 
dence and relative intensity of several superstructure peaks in the reciprocal lattice cell, we confirm 
antiferromagnetism with largely G-type character and T/v = 324 K, significantly below a structural 
phase transition at 375 K, contrary to previous reports. Evidence for a second transition, possibly a 
minority magnetic phase with C-type character is also reported with T/v = 260 K. The co-existence 
of the two magnetic phases in T-like BiFeOs and the difference in ordering temperatures between 
i?-like and T-like systems is explained through simple Fe-O-Fe bond distance considerations. 

PACS numbers: 75.25.-j,75.70.Ak,77.55.Nv,77.80.bn 



The oxide-perovskite BiFeOs is the only known room 
temperature multiferroic[3], with robust ferroelectricity 
below T c = 1103 K and long- wavelength cycloidal anti- 
ferromagnetism below T/v = 643 K in bulk form [4]. A 
flurry of recent work indicates that thin films can ex- 
hibit large ferroelectric polarization[5j[6], novel magneto- 
electric coupling effects [7 and a series of strain-induced 
phase transitions with the promise of compelling device 
applications [8-13 [ . Currently there is particularly strong 
interest in the epitaxially-stabilized tetragonal-like ('T- 
like') monoclinic phase of BiFeOs induced by large com- 
pressive strains. This Letter presents new experimental 
results on the nature of magnetism in essentially pure 
T-like films, and shows that the structural and antiferro- 
magnetic phase transitions are distinct from one another, 
contrary to previous reports [U|2]. 

Epitaxial strain is known to play a major role in deter- 
mining properties of thin-film multiferroic materials [T4]- 
[T6] . In BiFeOs, even the slight compressive strain af- 
forded by epitaxial growth on SrTiOs substrates changes 
the structure from bulk rhombohedral to 'itMike' mono- 
clinic with a - ratio near 1.03 [5]. The large room tem- 
perature ferroelectric moment of bulk BiFeOs is retained, 
and with increasing strain the crucial out-of-plane ferro- 
electric moment is enhanced [6 . Antiferromagnetic order 
is largely G-tvpe fT7lfT9] with a commensurability that 
depends on local structural considerations [T9] . 

When the compressive strain exceeds 4.5% a new phase 
appears with a - ratio near 1.25 [5J [9]. This 'T-like' phase 
can be thought of as a monoclinic distortion from a pre- 
dicted meta-stable PAmm tetragonal phase [20- 22 . It 



is now understood to be part of a sequence of phases 
in BiFeOs: with compressive strain the structural pro- 
gression is rhombohedral - (i?-like) monoclinic - (T-like) 
monoclinic - tetragonal , similar to what is observed 
at morphotropic phase boundaries in lead-based piezo- 
electric materials [13 . As has recently been pointed out 
by Damodaran et al. [23], this response to strain is very 
different from that of a Poisson-like deformation. In 
a Poisson-like elastic regime, coherent strain between 
a film a substrate is only maintained up to a critical 
thickness [24]. BiFeOs, in contrast, responds via polymor- 
phism similar to the epitaxial halides[25 or antimony, 
which can be epitaxially stabilized in the "grey tin" al- 
lotrope rather than the stable "white tin" phase in the 
bulk[26]. 

While the ferroelectric transition temperature, T c , of 
i?-like films decreases with compressive strain [28], its 
value is unknown for T-like BiFeOs. However, such films 
exhibit a structural phase transition between distinct T- 
like monoclinic phases at 375 K [TJ [TJ [30j [31] , leading to 
the anticipation of strong coupling effects as the various 
transition temperatures become comparable. 

Very little is known at this point about the magnetic 
properties of T-like phase of BiFeOs, in large part be- 
cause magnetometry studies of antiferromagnetic films 
are often inconclusive due to the strong diamagnetic 
response of substrate materials. First-principles cal- 
culations for films with 6% compressive strain found 
near-degenerate energy scales for G-type and C-type 
magnetic structures, where adjacent antiferromagnetic 
layers are stacked with spins antiparallel and parallel, 



respectively [22]. Such calculations [22| [32] predict that 
structures with a large ^ ratio, including the T-like 
phase, should exhibit ground states with C-type antifer- 
romagnetic order [32]. In contrast, a room temperature 
neutron diffraction measurement on one phase-pure T- 
like BiFeOs nmi revealed a magnetic peak that is consis- 
tent with a G-type structure [8]. Subsequent Mossbauer 
measurements have recently been used to infer an anti- 
ferromagnetic transition at T/v = 360 ± 20 KpQ, within 
error equal to the above-mentioned structural transition. 
This is remarkably lower than the transitions to G-type 
antiferromagnetism in i?-like films that occurs with T/v 
near the bulk value of 643 K [19] [28] . 

In this Letter, we report elastic neutron scattering 
measurements in epitaxially-stabilized T-like BiFeOs on 
LaA103 which go beyond previous studies in several im- 
portant ways and provide important new information 
about the nature of magnetism in this material. By 
measuring the temperature dependence of scattering at 
multiple half-integer Bragg positions, we confirm the 
existence G-type antiferromagnetism implied by earlier 
measurements [8] and provide the first estimate of Neel 
temperature with neutrons: T/v = 324 K. By addition- 
ally observing the T = 375 K structural transition in the 
same experiment, our data clearly demonstrate that the 
two transitions are in fact decoupled. We further present 
evidence for a second magnetic phase showing C-type an- 
tiferromagnetism at a lower temperature, demonstrating 
macroscopic phase co-existence in a T-like BiFeOs nmi 
that does not contain the .R-like polymorph. Through a 
comparison of scattering intensities at several magnetic 
Bragg peaks, it is inferred that spins in the G-type struc- 
ture lie in the plane of the film. Aided by Monte Carlo 
calculations, we show that the essential difference in the 
magnetism of i?-like and T-like films can be understood 
using straightforward considerations of Fe-O-Fe bond dis- 
tances. 

An epitaxial BiFeOs fil m of 300nm thickness was 
grown on (OOl)-oriented LaAlOs substrates via pulsed 
laser deposition. (Pseudo-cubic notations are used 
throughout this Letter.) Details about the growth tech- 
nique and conditions are described elsewhere [13], as are 
extensive x-ray scattering investigations of the quality 
and structural evolution of the film investigated in the 
current Letter pQ. 

The data in Ref. [1] (obtained on the same sample stud- 
ied here) show that the film consists only of T-like BiFeOs 
(^=1.23), a small (< 2%) amount of a novel polymorph 
that reversibly arises below the T =375 K transition due 
to tiling considerations, a negligible (<< 0.1%) of poly- 
crystalline bismuth oxide, and no measurable amount of 
the i?-like phase. This near phase purity contrasts some 
previous reports of BiFeOs films on LaAlOs, where a 
competing iiMike phase occupies up to 30-60% of the 
sample volume [TJ [9]. Nonetheless, it is in agreement 
with numerous reports in the literature of up to 99% 
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FIG. 1: Temperature dependence of the out-of-plane (dia- 
monds) and in-plane (circles) lattice parameters of the BiFeOs 
film. Closed symbols represent measurements taken with the 
high-temperature CCR apparatus upon warming from 40 K, 
and open symbols represent data taken with the furnace upon 
cooling from 500 K. A sharp change in the temperature de- 
pendence of c is seen at T* « 375 K (dashed line), without a 
concomitant change in a. This agrees with the T-like mono- 
clinic to T-like monoclinic structural phase transition recently 
observed with x-ray scattering[l] [T] [30] [31]. Solid lines are 
guides to the eye. 



phase pure T-like BiFeOs in samples of thickness lOOnm 
or greater [HE] [30]. 

Neutron scattering measurements were performed us- 
ing the HBla Triple- Axis Spectrometer at the High-Flux 
Isotope Reactor in an elastic configuration with a fixed 
incident energy T^=14.7meV, pyrolytic graphite (PG) 
(002) monochromator and analyzer, instrument collima- 
tions 48 / -48 / -40 / -60 / , and higher order contamination re- 
moved by standard PG filters. The film was aligned in 
the (H H L) scattering plane. For measurements at tem- 
peratures ranging from 10 K to 375 K, the sample was 
mounted in a closed cycle refrigerator (CCR) in a can 
with thermal contact assured by He exchange gas. Addi- 
tional measurements from 40 K to 435 K were done using 
a high-temperature CCR which lacked exchange gas and 
from 300 K to 500 K using a furnace. Where presented 
below, data from the furnace were scaled to match the 
data from the CCR as laid out in the supplementary ma- 
terial. 

Neutron measurements confirmed that the BiFeOs fil m 
had a nearly tetragonal structure with the in-plane lat- 
tice constants epitaxially matched to the substrate. The 
room temperature lattice constants were a = 3. 79 A, c 
= 4.66A- consistent with the - ratio for the T-like phase 
reported by Christen et al [13] and results from x-ray scat- 



tering on this sample [1. Measurements of the film peak 
at (001) using both the high-temperature CCR and the 
furnace were used to track the evolution of the out-of- 
plane lattice parameter, c, with temperature. The re- 
sults, shown in Fig. [T] reveal a significant evolution of c 
over the entire temperature range investigated, while the 
in-plane lattice parameter remained perfectly epitaxially 
matched to the substrate within resolution. A sharp cusp 
in the c-temperature curve at T* « 380 K reflects the ex- 
istence of a structural phase transition between different 
T-like phases of BiFeOs, as reported in recent x-ray and 
Raman scattering studies [TJ [TJ [30j [31] . 

At low temperatures, neutron scattering showed super- 
lattice peaks in the film at half-integer Bragg positions 
for the majority structural phase. The large ^ ratio of the 
T-like phase allows us to clearly distinguish between sig- 
nal from it and one from, for example, a possible i?-like 
impurity phase that might have escaped our attention 
in x-ray scattering measurements. Most of the observed 
peaks were characteristic of G-type order: ( \\ \ ) , ( \ \\ ) , 
(~|), (||^) and (|||)- In addition, a peak was ob- 
served at the C-type wave- vector (^|0). No significant 
signal was apparent at (00^) or (11 1) positions, elimi- 
nating A-type antiferromagnetism. Fig. [2] shows radial 
scans across four peaks at temperatures above (370 K) 
and below (37 K) the Neel temperatures determined be- 
low. All panels show a significant increase of scattering 
intensity at low temperature. The weak peaks remaining 
at 370 K in panels (a) and (c) might arise from resid- 
ual magnetic short range order or weak structural distor- 
tions. The temperature independent background peaks 
in panels (b) and (c) are powder diffraction lines of un- 
known origin; it has been confirmed that the background 
peak in panel (d) is dependent on the choice of sample 
environment. 

Temperature dependent x-ray scattering confirmed 
that there is no structural transition in this sample be- 
tween 300 K and 375 K[l . The neutron data presented 
in Fig. [I] suggest no major modifications of the film struc- 
ture over the entire temperature range investigated. At 
low temperatures, the observed intensity at (|||) is 
slightly less than 10% of that at (\\\), consistent with 
expectations based on the magnetic form factor of Fe 3+ . 
Moreover, the intensity ratios of the half-integer peaks 
are inconsistent with any of the known common struc- 
tural distortions observed in cubic perovskites. The lim- 
ited data set available does not permit a complete re- 
finement of the magnetic structure, however the relative 
intensities of different {\\L) peaks imply that moments 
of the G-type structure are confined to the plane of the 
film. To within the resolution of the present measure- 
ments, no modulation of the G-type structure is visible. 
As mentioned above, the emergence of a peak at the (\\ 
0) Bragg position is at odds with a simple G-type pic- 
ture. We further point out that the scattering at the 
(t^0) position is much too intense to be associated with 
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FIG. 2: Elastic neutron scattering data at four half-integer 
film peaks at T = 370 K (open circles) and 37 K (closed 
diamonds) measured using the high-temperature CCR envi- 
ronment. Clear increased intensity at low T is seen at peaks 
associated with G-type antiferromagnetic order. Shown here 
are: (a)(|||), (b)(||§) and (c)(|||). Additional intensity 
at low T is also observed at (d) (||0), which may indicate 
a small co-existing volume of C-type antiferromagnetism (see 
discussion in text). Arrows in panels (b) and (d) denote po- 
sitions of temperature-independent background peaks. Solid 
lines are the results of fitting to simple Gaussians. 



the minority structural phase identified by x-ray scatter- 
ing, which would appear an order-of-magnitude weaker in 
neutron scattering experiments. Since only one peak has 
been detected, one cannot definitively rule out a struc- 
tural phase transition, albeit with no visible effect in the 
lattice parameters in Fig. [I] More likely is the onset of 
a co-existence of two distinct magnetic ordering patterns 
at lower temperatures, both within the majority T-like 
structural phase. This latter possibility is consistent with 
the near-degeneracy of G-type and C-type structures in 
the ab-initio calculations [22]. A comparison of the rela- 
tive intensities of (\\\) and (||0) peaks suggests that 
15-30% of the sample volume contains C-type antiferro- 
magnetism in this scenario, depending on the orientation 
of the ordered spins. As the widths of the superlattice 
peaks measured in low-order Brillouin zones are compa- 
rable to the width of the substrate (2 2 0) peak, it is not 
possible in the present experiment to infer further infor- 
mation about the magnetic domain sizes or correlation 
lengths. 

The magnetic transition temperature of the G-type 
order was determined from the measured intensity of 
the scattering at (|||), as shown in the main panel of 
Fig. [3] Data from both the CCR and furnace were fit phe- 
nomenologically to a power-law temperature dependence 
in the range 10K<T< T/v. Particular attention was also 
given to the fits of background above the magnetic transi- 
tion, as a small amount of temperature dependence was 
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FIG. 3: Plot of the neutron scattering intensity at the (1/2 
1/2 1/2) peak position as a function of temperature, measured 
with a furnace (closed diamonds) and a low-temperature CCR 
(open circles) upon cooling. The data from the furnace were 
adjusted by scaling using the ratio of measured Bragg peak 
heights and adjusting the background to match the range of 
data overlapped with the CCR. The solid (dashed) line is a fit 
assuming sloped (flat) background. Inset shows the evolution 
of the fitted (1/2 1/2 0) peak height over background scaled 
to match counting time and conditions of main panel. 



observed up to 500 K, possibly associated with critical 
correlations above a continuous phase transition. Specif- 
ically, fits were performed with either linearly sloping or 
flat background. The former fit, shown as a solid line 
in Fig. |3j provides the superior description of the data 
and implies a transition temperature of T/v = 323 ± 4 
K. We believe this to be the most reliable estimate of 
T/v, however it should be noted that the fit with back- 
ground held constant (dashed line in Fig. [3]) yields T/v 
= 341 ± 3K. The difference in these values reflects the 
magnitude of systematic error that might arise from dif- 
ferent choices of background, and in fact these values may 
be taken as upper and lower bounds of reasonable esti- 
mates of T^. For example, if a flat background is used 
and temperatures where the curve is concave-up are ex- 
cluded, fits imply T/v = 327 ±5K. All values represent a 
significant decrease from the Neel temperatures reported 
for bulk BiFe0 3 (643 K)[3] and iMike films (600 - 650 
K)[28|. The magnetic transition is also clearly distinct 
from the structural transition identified in Fig. [T] in con- 
trast to some previous reports [TJ [2]. The temperature 
dependence of the (|^0) peak intensity, shown in the in- 
set of Fig. [3j suggests that the C-type order is seen only 
at temperatures below Tn,c — 264 ± 12 K. (Error bars 
are statistical and do not include possible sources of sys- 
tematic errors.) No change is seen in the intensity of the 
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FIG. 4: The convergence tendency of the ordering tempera- 
tures T/v, deduced from Monte Carlo calculations of staggered 
magnetization, susceptibility and heat capacity in /Mike and 
T-like BiFeOs films. Details about the model used in the 
calculations are given in the main text. The star represents 
the value of T/v determined with neutrons. Also depicted are 
the magnetic ordering patterns observed in itMike and T-like 
BiFeOs, with ordered moments chosen to lie in the a-b plane. 



(\\\) scattering at T N ^ C - 

Simple considerations can provide insight into the dif- 
ferent magnetic properties of the iiMike and T-like phases 
of BiFeOs. Transition metal pseudopotential theory [33 
predicts that the hopping matrix elements t p d between 
the d-orbitals of Fe and the p-orbitals of O scale with 
the Fe-0 interatomic distance r as -^ . Since perturba- 
tion theory near the atomic limit predicts that the Fe-Fe 



Heisenberg superexchange J depends on t p d via J ex t 4 



pdi 



then overall J varies as J oc r . Using the known 
structures [13J and normalizing by the in-plane exchange 
constant of .R-like BiFeOs, this criterion implies that the 
i?-like superexchange is characterized by the parameters 
J<E=«/y=1.0 and out-of-plane J 2 =0.87, while the T-like 
superexchange is characterized by J x =1.43, J y =2Al 1 and 
J z =0.10. Although the precise results will be modified 
by additionally considering true bond angles, it is imme- 
diately apparent that the interactions in T-like BiFeOs 
are expected to be much more two-dimensional than the 
interactions in the i?-like material. Even though the in- 
plane values of exchange are higher for the T-like phase, 
our data imply that the net result is a suppression of 
magnetic order to much lower temperatures. 

To explore this idea further, we have performed Monte 
Carlo simulations of the two cases utilizing an anisotropic 
classical Heisenberg Hamiltonian with exchange cou- 
plings in the ratios discussed above. Calculations were 
carried out on a set of LxLxL sized tetragonal lattices, 



with transition temperatures determined by examining 
the staggered magnetization, susceptibility and heat ca- 
pacity. The precise magnitude of the superexchange was 
scaled to match the observed Neel temperature of the 
i?-like BiFe0 3 films, 643 K. Fig. [4] plots the variation of 
T/v for the two cases with the Monte Carlo system size, 
L. The convergence to the bulk limit is much slower for 
the T-like system, as expected because of its lower effec- 
tive dimensionality. However, the present Monte Carlo 
results are sufficient to show that the ordering for the T- 
like phase will occur at a substantially lower temperature 
than for the i?-like phase, in qualitative agreement with 
the experimental observations. 

The weak inter-planar coupling of T-like BiFeOs also 
offers some insight into the possible co-existence of G- 
type and ( \\ 0) C-type antiferromagnetism at low tem- 
peratures. Both structures (illustrated in Fig. [4| are 
characterized by an antiferromagnetic arrangement of 
spins in the a-b plane, where the three-dimensional G- 
type structure is realized by antiferromagnetic stacking 
and C-type by ferromagnetic stacking along (00L). With 
very weak inter-planar interactions, one can easily sup- 
pose that small differences in local structure or impurities 
can lead to either state, and this is consistent with the 
findings of ab-initio calculations that the energy differ- 
ences are very small and may in fact favor C-type order. 

In summary, the results presented here shed light on 
the physics of the epitaxially-stabilized T-like monoclinic 
phase of BiFeOs which only exists in samples with high 
levels of compressive strain. Our film grown on LaAlOs 
is nearly crystallographically phase pure, but gives evi- 
dence for phase co-existence of G-type and C-type anti- 
ferromagnetic order parameters. We estimate magnetic 
ordering temperature T/v = 323 ± 4 K for the G-type 
order, which is distinct from the structural phase transi- 
tion seen at higher temperature and more than a factor 
of two lower than similar order seen in i?-like films. We 
attribute this large reduction mainly to a reduced dimen- 
sionality associated with increased Fe-Fe distances in the 
out-of-plane direction. 
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SUPPORTING INFORMATION 



SAMPLE ENVIRONMENTS AND THE SCALING 
OF DATA SETS. 



SAMPLE CHARACTERIZATION 

The epitaxial BiFeOs fil m of the current study was 
grown on (OOl)-oriented LaAlOs substrates via pulsed 
laser deposition. Thickness was determined to be 300nm 
(± 20%) from the spacing of x-ray thickness fringes and 
growth rate calculations. The sample was then charac- 
terized using a PANalytical X'Pert 4-circle thin-film x- 
ray diffractometer using Cu-Ka radiation and with the 
sample mounted on an Anton-Paar heated stage. The 
resultant data have largely been published and discussed 
at length in Reference 1. A single x-ray pattern taken 
at room temperature is reproduced in Fig. [I] for the pur- 
pose of demonstrating sample quality. As is easily seen, 
x-ray scattering data indicate that the film consists al- 
most entirely of the T-like monoclinic phase of BiFeOs, 
with ^=1.23. Impurity phases are reduced by nearly two 
orders of magnitude (note the logarithmic scale in Fig. [I]), 
and consist of a small (< 2%) amount of a novel T-like 
polymorph of BiFeOs that reversibly arises below the T 
=375 K transition from tiling considerations (see Ref. [T]) 
and a negligible (<< 0.1%) amount of poly crystalline 
bismuth oxide. There is no measurable amount of the 
i?-like phase of BiFeOs, which would manifest as peaks 
nearer to those of LaAlOs with -=1.03. 
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FIG. 1: X-ray diffraction pattern from the sample in the cur- 
rent study, plotted on a logarithmic scale. Peaks arising from 
the film (BFO), the substrate (LAO) and the small BiO impu- 
rity have been labeled explicitly. Also present are peaks from 
the sample holder and a new T-like polymorph, labeled with 
black dots and stars, respectively. Data have been published 
previously and discussed at length in Ref. 1 



Neutron scattering data presented in the main article 
were obtained using three separate sample environments 
for controlling temperature over the range 10-500K: a 
high-temperature closed-cycle refrigerator (CCR) over 
the range 40-400 K, a high-temperature furnace over the 
range 300-500 K, and a low-temperature CCR with ex- 
change gas which covered the range 10-370 K. It was 
demonstrated in Figure 1 of the main article that both 
data from the furnace and high-temperature CCR re- 
flect the existence of the structural phase transition pre- 
viously observed by x-ray scattering and can be used to 
independently verify the structural transition tempera- 
ture of T s ~380 K. As the measurements of the c-axis 
lattice parameter were measured on warming in the high- 
temperature CCR and upon cooling in the furnace, this 
further demonstrates that there is minimal hysteresis for 
this transition in our sample. In neither of these high- 
temperature sample environments is there any signature 
of magnetism at or appearing immediately below the 
structural transition temperature. 

A measure of the Neel temperature for the G-type an- 
tiferromagnetically ordered phase can be determined in- 
dependently using each of the three data sets and the 
estimates are found to be consistent with one another. 
One example is shown in Fig. [2|a), where plotted ver- 
sus temperature is the integrated intensity of the (|||) 
magnetic Bragg peak as measured via radial scans on 
warming in the high-temperature CCR. It is clear from 
this data that T/v is above room temperature and no 
larger than 340 K, consistent with the analysis presented 
in the main text (see Fig. 3 of the main Letter). To check 
for possible hysteresis effects, which would imply a lower 
T/v, we plot in Fig. |2jb) the peak height extracted from 
the same radial scans and compare to measurements of 
the scattering intensity measured while cooling in the 
same sample environment. The two data sets agree over 
an extended temperature range, implying that for this 
transition as well hysteresis is minimal. 

For direct comparison of the different data sets, one 
can also scale data according to: 



I (T) = a(I(T) - Bd 



instrument 



) + /?, 



(i) 



where a accounts for differing effective scattering vol- 
ume between separate measurements and BGi nstr ument 
is associated with sample-independent sources of scatter- 
ing background. Fig. [3] demonstrates the result scaling 
data from all three sample environments to match those 
taken in the low-temperature CCR. Here, BGi nstr ument 
was chosen to be the temperature independent part of 
the Q-independent background extracted from radial 
scans across the magnetic Bragg peak; this is a valid as- 
sumption since incoherent scattering in both BiFeOs and 
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FIG. 2: (a) Integrated intensity of the (| 1 1) magnetic Bragg 
peak, as measured via radial scans taken on warming using the 
high-temperature CCR sample environment, (b) Peak height 
(red stars) and the Q-independent background (black circles) 
extracted from the same scans and plotted as a function of 
temperature. Also plotted is the neutron scattering intensity 
at the (|||) taken in the same sample environment while 
cooling from 400 K, demonstrating the minimal amount of 
hysteresis observed in the current sample. 



LaAlOs is minimal. The scaling factors a were chosen by 
comparing nuclear Bragg peak intensities of the film near 
room temperature, and /3 was chosen such that magnetic 
data were consistent at T=370K (T N < T < T* tructural ). 
Solid and dashed lines are the fits described in the main 
text. Of immediate note is that at high temperature all 
of the scattering data are consistent and imply a similar 
magnetic transition temperature. At lower temperatures, 
the scattering measured using the high-temperature CCR 
is systematically lower than what is measured in the low- 
temperature CCR. We attribute this to inefficient cooling 
of the sample in the high-temperature displex, expected 
in this sample environment at low temperatures due to 
the lack of exchange gas. 



DETERMINING THE NEEL TEMPERATURE 

When estimating T/v in the main article, our analy- 
sis focussed on data obtained using the low-temperature 
CCR instead of the high-temperature CCR data plotted 
in Figure 2. This was because the low-temperature CCR 
data was deemed more reliable, due to both the large 
uncertainty in temperature during the high-T CCR mea- 
surements and the existence of helium exchange gas in 
the low-temperature apparatus. We additionally note the 
temperature-dependent scattering at the (|||) position 
apparent in all of our data up to 500 K. This might rep- 
resent short-range critical correlations, which must exist 
above a continuous phase transition. In modeling the 
temperature dependence of this short-range scattering, 
we have considered the cases of a linearly varying back- 
ground, a temperature independent (flat) background 
and a flat background with the critical region excluded 
from fits. 

In this section, we demonstrate directly that neither 
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FIG. 3: Plot of neutron scattering intensity at the (|||) 
position of the BiFe03 him as a function of temperature, 
including data taken in the low-T CCR (circles), the high- 
temperature furnace (diamonds), all three sample environ- 
ments employed in this experiment. All data were scaled to 
that in the low-temperature CCR, as laid out in the Supple- 
mentary text. 



the choice of data set nor specific modeling of the back- 
ground temperature dependence has a significant effect 
on conclusions of the current paper. Data from differ- 
ent sample environments were fit phenomenologically, as- 
suming the follow form for magnetic scattering: 



I(T) = Jo 



LN 



T N 



2/3 



■ Ibg, 



(2) 



with Iq and j3 as fit parameters and Ibg determined 
by data at temperatures T>380 K. Table 1 displays the 
results of this analysis, giving the inferred value of T/v 
and quality of fit parameter y 2 for each scenario. Error 
bars reported here are statistical, arising from fits and 
not accounting for possible systematic effects. This is an 
especially relevant concern when considering data taken 
in the high-temperature CCR, as the measurements were 
almost certainly taken in non-equilibrium conditions due 
to equipment failure. One can see, however, that each fit 
in Table 1 implies a transition temperature in the range 
322 < T/v < 356 K. Even the highest estimate of T/v is 
found to be 5 standard deviations below the known struc- 
tural transition temperature at 375 K. If one excludes the 
less reliable high-T CCR data, this range is even more 
restrictive: 322< T/v < 341 K. For every choice of data 
subset, the fits assuming sloping background were found 
to describe the data best. It is worth further noting that 
for these fits, there is very little variation in inferred T/v 



among data sets. Our best estimate of Neel temperature 
is thus determined to be T/v= 323±4 K. 



Data set 


BG function 


T N (K) 


x 2 
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325+3 


2.14 


L+F 
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324+4 


1.90 
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s 


323+7 


0.99 


H+F 


s 


322+5 


0.87 


H+L+F 


s 


324+3 


1.77 


L 


fl 


342+3 


2.78 


L+F 


fl 


341+3 


2.61 


H 


fl 


356+4 


1.23 


H+F 


fl 


356+4 


1.36 


H+L+F 


fl 


344+2 


2.4 


L 


fl-x 


327+5 


1.67 


L+F 


fl-x 


327+5 


1.69 


H 


fl-x 


329+12 


1.20 


H+F 


fl-x 


329+12 


1.29 


H+L+F 


fl-x 


325+4 


1.70 



TABLE I: The results of fitting scattering intensity at the 
(|||) versus temperature, using data from different sam- 
ple environments and under different assumptions about the 
high-temperature background. Considered are fits to the 
low-temperature CCR alone (L), the L-CCR plus the fur- 
nace(L+F), the high-temperature CCR alone (H), the high- 
temperature CCR plus the furnace (H+F) and the entirety of 
available data (H+L+F). Backgrounds are considered to be 
either sloping (s) with temperature, flat (fl) or flat with the 
critical temperature region excluded from the fit (fl-x). 



Electronic address: christenhm@ornl.gov 
[1] W. Siemons, M. D. Biegalski, J. H. Nam, and H. M. Chris- 
ten, Appl. Phys. Express 4, 095801 (2011). 



